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SECTION 1 INTRODUCTION 

O'Keefe (1959) tested the basic hypothesis of geodesy by 
using the harmonics of the earth ,?s gravity field from satellite 
measurements, He pointed out that these harmonics may be 
supported by hydrodynamic forces resulting from mantle convection. 
Runcorn (1964, 1967) developed a theory connecting mantle con- 
vection and the geoid. 

Although numerical calculations of thermal convection models 
of the mantle have been carried out by Turcotte et al. (1973) and 

a. 

by McKenzie et al. (1974) , the non-linear nature of mantle 
convection and the wide range of viscosities present formidable 
obstacles to the thorough knowledge of convection patterns. It 
may well be that unless simple models of mantle flow are first 
understood, complex cases would not be interpretable in terms 
of the contribution to the dynamics that each case provides. 
Therefore, it seems necessary, and perhaps even desirable, to 
apply the relatively simple model of mantle convection developed 
by Runcorn (1964, 1967). In this paper, we propose to detect 
and determine the tectonic forces or stresses due to the small- 
scale mantle flow under the South American plate by utilizing the 
harmonics of the geopotential field model. 
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SECTION 2 THEORY OF MANTLE CONVECTION AND THE GEOID 


Given a set of fully normalized harmonic coefficients, 

C ^ and . determined from satellite measurements and surface 
n,m n,ra 

gravity observations, the relationship between the mantle 
convection pattern at the upper boundary and the geoid may be 
described by (Runcorn, 1964, 1967): 
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where 

M = mass of the earth, 

y * coefficient of viscosity, 

a * radius of outer spherical surface of the flowing region, 

(cos?) = associated Legendre polynomial of degree n, 
order m and argument cos , 

X = the longitude, 

c = n - c 1 

^ 2,0 '“ 2 , 0 c 2 , 0 ' 

S = S 
n,m n,m 

C 2 , 0 and 

C^ q = coefficients characterizing the ellipsoid of revolution, 

W = a poloidal function of convection currents 

g = gravitational acceleration, 

a' = radius of the earth, 

? = the co-latitude. 
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SECTION 2 (Continued') 


5 4,0 * C 4 , 0 “ C 4 , 0 ' 

C - * C „ for nj*2 or 4 and m«0. 

n,m n,m 

In th® derivation of eq. 1, a Newtonian viscosity is attributed 
to the mantle and laminar viscous flow is assumed. 
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SECTION 3 SUBCRUSTAL STRESSES EXERTED BY 
CONVECTION CURRENTS 


If the crust above the flowing mantle is assumed to be 

an elastic shell, the tangential component of the velocity 

vanishes at r»a, i.e., (dW /dr) I . ■ 0. Then the stress 

n r*a 

components exerted by the convection flow on the crust in the 
eastward and northward directions are determined by: 
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SECTION 4 DATA SETS 


The harmonic coefficients C _ and S_ _ of the geopotential 

n,m n,m 

field model in eqs. 2 and 3 have been significantly improved in 
recent years, principally because of the improved satellite 
observations and surface gravity measurements (Smith et. al. , 1976) . 
The low-degree harmonics for n < 12 may reflect a large-scale 
mantle flow system (Runcorn 1967; Richter and Parsons, 1975). 

In this paper, the high-degree harmonics, n > 13, are assumed to 
describe the small-scale mantle convection patterns. The GEM-10 
gravity model used here is based on the work by Lerch et. al. (1977) 
The input data used in the derivation of this model is made up of 
840,000 optical, electronic and laser observations and 1,656 
5° x 5° mean free-air anomalies. Although there remain some 
statistically questionable aspects of the high-degree harmonics, 
it seems appropriate now to explore their implications for the 
tectonic forces or stress field under the crust. 
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SECTION 5 CALCULATIONS AND RESULTS 


0 


In order to compute the magnitudes and directions of 
stresses exerted by the convection currents on the crust , the 
following recursive relationships among the Legendre and associated 
Legendre functions were employed (Arfken, 1966) : 

(2n+l) cos?P® (cos?) * (n+1) pJ? , . (cose) + nP® .(cos?) (4) 

n n+i n-i 


S (cos?) - p“ + }(cos?) - (2n+l) sin? P® (cos?) 

n+i n-x 
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= (2n+l) sin? g-y- P“ l (cos?) + (2n+l)cos? P™ (cos?) 
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SECTION 5 (Continued) 


As each function was computed, it was normalized using 


(Mather, 1971): 



where 


(2-6 Qm ) (2n+l) (n-rn) ! 


(n+m) ! 


r 


p JJ(u) 



jo if m^O 
( 1 if m«0 


Due to the slow computer job turn-around at the GSFC 360/95, 
the program SCP that computes the stress components as a function 
of latitude and longitude was coded so that each job computed the 
stress components for a small region and the results from each 
job were later merged. This merged result was then used as input 
to the stress plotting program SPP which generated the necessary 
data for producing a picture display on the SD 4060. 

The figure shows the stress pattern corresponding to the 
mantle flow system under the South American Plate. As shown in 
the figure, in the vicinity of Bolivia and off the coast of 
Bolivia, these regions are under tension. These results are in 
excellent agreement with the geoidal highs that have been concluded 
by Marsh and Chang (1979). 
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